The clustering instability of inertial particles spatial distribution in turbulent flows 
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A theory of clustering of inertial particles advected by a turbulent velocity field caused by an 
instability of their spatial distribution is suggested. The reason for the clustering instability is a 
combined effect of the particles inertia and a finite correlation time of the velocity field. The crucial 
parameter for the clustering instability is a size of the particles. The critical size is estimated for 
a strong clustering (with a finite fraction of particles in clusters) associated with the growth of 
the mean absolute value of the particles number density and for a weak clustering associated with 
the growth of the second and higher moments. A new concept of compressibility of the turbulent 
diffusion tensor caused by a finite correlation time of an incompressible velocity field is introduced. 
In this model of the velocity field, the field of Lagrangian trajectories is not divergence-free. A 
mechanism of saturation of the clustering instability associated with the particles collisions in the 
clusters is suggested. Applications of the analyzed effects to the dynamics of droplets in the turbulent 
atmosphere are discussed. An estimated nonlinear level of the saturation of the droplets number 
density in clouds exceeds by the orders of magnitude their mean number density. The critical size 
of cloud droplets required for clusters formation is more than 20/im. 

PACS numbers: 47.27.Qb, 05.40.-a 



I. INTRODUCTION 

Formation and evolution of aerosols and droplets in- 
homogeneities (clusters) are of fundamental significance 
in many areas of environmental sciences, physics of the 
atmosphere and meteorology (e.g., smog and fog for- 
mation, rain formation), transport and mixing in in- 
dustrial turbulent flows (like spray drying, pulverized- 
coal-fircd furnaces, cyclone dust separation, abrasive 
water-jet cutting) and in turbulent combustion (see, e.g., 
§, |, |, 1, |, |,§ 1). The reason is that the direct, hy- 
drodynamic, diffusional and thermal interactions of par- 
ticles in dense clusters strongly affect the character of the 
involved phenomena. Thus, e.g., enhanced binary colli- 
sions between cloud droplets in dense clusters can cause 
fast broadening of droplet size spectrum and rain for- 
mation (see, e.g., ||). Another example is combustion of 
pulverized coal or sprays whereby reaction rate of a single 
particle or a droplet differs considerably from a reaction 
rate of a coal particle or a droplet in a cluster (see, e.g., 

Analysis of experimental data shows that spatial distri- 
butions of droplets in clouds are strongly inhomogeneous 
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(see, e.g., ]ll], [l2[ [l3| 0). Small-scale inhomogeneities 
in particles distribution were observed also in laboratory 
turbulent flows (r| ^ |§. 

It is well-known that turbulence results in a relaxation 
of inhomogeneities of concentration due to turbulent dif- 
fusion, whereas the opposite process, e.g., a preferen- 
tial concentration (clustering) of droplets and particles 
in turbulent fluid flow still remains poorly understood. 

In this paper we suggest a theory of clustering of parti- 
cles and droplets in turbulent flows. The clusters of par- 
ticles are formed due to an instability of their spatial dis- 
tribution suggested in Ref . [[l9| and caused by a combined 
effect of a particle inertia and a finite velocity correlation 
time. Particles inside turbulent eddies are carried out 
to the boundary regions between them by inertial forces. 
This mechanism of the preferential concentration acts in 
all scales of turbulence, increasing toward small scales. 
An opposite process, a relaxation of clusters is caused by 
a scale-dependent turbulent diffusion. The turbulent dif- 
fusion decreases towards to smaller scales. Therefore, the 
clustering instability dominates in the Kolmogorov inner 
scale r], which separates inertial and viscous scales. Expo- 
nential growth of the number of particles in the clusters 
is saturated by their collisions. 

In our previous study Jl!| we suggested and analyzed 
qualitatively an idea that inertia of particles may lead to 
their clustering. Later this idea was questioned by our 
quantitative analysis [^0| ^lj of the Kraichnan model of 
turbulent advection of particles by the delta-correlated in 
time random velocity field. It was proved that the clus- 
tering of inertial particles does not occur in the Kraich- 
nan model. The latter result may be considered as coun- 
terexample. 
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The main quantitative result of the theory of clustering 
instability of inertial particles, suggested in this paper, 
is the existence of this instability under some conditions 
that we determined. We showed that the inertia of the 
particles is only one of the necessary conditions for par- 
ticles clustering in turbulent flow. In the present study 
we found a second necessary condition for the clustering 
instability: a finite correlation time of the fluid velocity 
field which in the suggested theory results in a nonzero di- 
vergence of the field of Lagrangian trajectories. This time 
is equal zero in the above mentioned Kraichnan model 
(sec [Q) which was the reason for the disappearance of 
the instability in this particular model. 

In this study we used a model of the turbulent veloc- 
ity field with a finite correlation time which drastically 
changes the dynamics of inertial particles. In the frame- 
work of this model of the velocity field we rigorously de- 
rived the sufficient conditions for the clustering instabil- 
ity. We demonstrated the existence of the new phenom- 
ena of strong and weak clustering of inertial particles in 
a turbulent flow. These two types of the clustering in- 
stabilities have different physical meaning and different 
physical consequences in various phenomena. We com- 
puted also the instability thresholds which are different 
for the strong and weak clustering instabilities. 



II. QUALITATIVE ANALYSIS OF STRONG 
AND WEAK CLUSTERING 

A. Basic equations in the continuous media 
approximation 

In this study we used the equation for the number den- 
sity n(t, r) of particles advected by a turbulent velocity 
field u(t, r): 



dn(t, r) 
dt 



V • [n(t, r) v(t, r)] = D An(t, r) , (2.1) 



where D = kT /Qitvpa is the coefficient of molecular 
(Brownian) diffusion, v is the fluid kinematic viscosity, 
p and T are the fluid density and temperature, respec- 
tively, a is the radius of a particle and k is the Boltz- 
mann constant. Due to inertia of particles their velocity 
v(t,r) ^ u(t,r), e.g., the field v(t,r) is not divergence- 
free even for div u = (see |l^]). Equation ( |2.1[ ) implies 
conservation of the total number of particles in a closed 
volume. Consider 



Q(t,r) = n(t,r) 



(2.2) 



the deviation of n(t, r) from the uniform mean number 
dens ity o f particles n. Equation for Q(t, r) follows from 
Eq. Q: 



dO{t, r) 
dt~ 



[v(t,r)-V]@(t,r) (2.3) 
-Q(t, r) divv(t, r) + DAQ(t, r) . 



Here we assumed that the mean particles velocity is zero. 
We also neglected the term oc n div v describing an effect 
of an external source of fluctuations. This term does not 
affect the growth rate of the instability. In the present 
study we investigate only the effect of self-excitation of 
the clustering instability, and we do not consider an ef- 
fect of the source term on the dynamics of fluctuations. 
The source term cx n div v causes another type of fluctu- 
ations of particle number density which are not related 
with an instability and are localized in the maximum 
scale of turbulent motions. A mechanism of these fluctu- 
ations is related with perturbations of the mean number 
density of particles by a random divergent velocity field. 
The magnitude of these fluctuations is much lower than 
that of fluctuations which are caused by the clustering 
instability. 

In ou r qualitative analysis of the problem we use 
Eq. fl2.3| ) written in a co-moving with a cluster reference 
frame. Formally, this may be done using the Bclinicher- 
L'vov (BL) representation (for details, see |2^, Q). Let 
£~L{to, r \ t) will be Lagrangian trajectory of the reference 
point [in the particle velocity field v(t, r)} located at r at 
time to and pi,(to,r\t) be an increment of the trajectory: 

PL(t ,r\t) = f v[T,&(to,r\r)]dT, (2.4) 

Jt 

£ h (t ,r\t) = r + p L (t ,r\t) . 

By definition pl (to , f \ to ) =0 and £ (to , r | to) = r . Define 
as ro a position of a center of a cluster at the "initial" 
time to — (for the brevity of notations hereafter we 
skip the label to ) and consider a "co-moving" reference 
frame with the position of the origin at Co(t) = C( r o| 
Then BL velocity field v(ro\t,r) and BL velocity differ- 
ence W(ro\t,r) are defined as 

v(r \t,r) = v[t,r + p h (r \t)], (2.5) 
W(r \t,r) = v(r \t,r)-v(ro\t,ro) ■ (2.6) 

Actually the BL representation is very similar to the La- 
grangian description of the velocity field. The difference 
between the two representations is that in the Lagrangian 
representation one follows the trajectory of every fluid 
particle r + p-^(r,to\t) (located at r at time t = to), 
whereas in the BL-representation there is a special ini- 
tial point ro (in our case the initial position of the center 
of the cluster) whose trajectory determines the new co- 
ordinate system (see |^3|, |24|). With time the BL-field 
v(ro\t, r) becomes very different from the Lagrangian 
velocity field. It must be noted that the simultaneous 
correlators of both, the Lagrangian and the BL-velocity 
fields, are identical to the simultaneous correlators of the 
Eulcrian velocity v(r, t). The reason is that for station- 
ary statistics the simultaneous correlators do not depend 
on t, and in parti cula r one can assume t = to- 

Similar to Eq. (2.5) let us introduce BL representation 
for9(£,r): 



O(r \t,r) = e[t,r + p h (r \ 



(2.7) 
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In BL variables denned by Eqs. Eq. ( gj ) 

reads: 



dO(r \t,r) 
dt 



[W r (r o |t,r)-V]0(ro|t,r) (2.8) 



-0(r o |t, r) div W(r \t, r) + D A <d(r \t, r) 



The diff erence between Eqs. (2.3) and (2.8) is that 



Eq. (2.8) involves only velocity difference ( |2.6[ ) in which 



the velocity v(ro\t, t*q) of the cluster center is subtracted. 



B. Rigid-cluster Approximation 

Consider qualitatively a time evolution of different sta- 
tistical moments 

of the deviation Q(t, r) defined by 



M 9 (t) = (m,r)\*) 



(2.9) 



assuming that at the initial time, t = 0, the spatial distri- 
bution of particles is almost homogeneous, all moments 
A4 q (0) are small, where (•)„ denotes the ensemble aver- 
aging over random velocity field v. In order to eliminate 
the kinematic effect of sweeping of the cluster as a whol e 



we consider Eq. fl2.3| ) in the BL-representation, Eq. (|2 



Since the simultaneous moments of any field variables in 
the Eulerian and in the BL-representations coincide, the 
moments M q (t) can be written as 



M„(t) = <|e(r |t,r)|'> 



(2.10) 



Our conjecture is that on a qualitative level we can con- 
sider the role of each term in the Eq. ( |2.8| ) separately, as- 
suming some reasonable, time independent, frozen shape 
9(x) of a distribution Q(ro\t, r) inside a cluster: 



e(r |t, 



A(t)6 



(2.11) 



Here A(t) is time-dependent amplitude of a cluster, t c \ 
is the characteristic width of the cluster. Shape func- 
tion 0(x) may be chosen with a maximum equal one at 
x = and unit width. Real shapes of various clusters 
in the turbulent ensemble are determined by a compe- 
tition of different terms in the evolution equation (2 



However, we believe that particular shapes affect only 
numerical factors in the expression for the growth rate 
of clusters and do not effect their functional dependence 
on the parameters of the problem which is considered in 
this subsection. 



Effect of turbulent diffusion 



The advective term in the LHS of Eq. (2.8) results in 



turbulent diffusion inside the cluster. This effect may be 
modelled by renormalization of the molecular diffusion 
coefficient D in the right hand side (RHS) of Eq. fl2.8j) 



by the effective turbulent diffusion coefficient Dt with a 
usual estimate of Dt' 



D 



D + Dt, D t ~ tdVd/3 



(2.12) 



Hereafter v c \ is the mean square vel ocit y of particles at 
the scale £ c \. Instead of the full Eq. (2.8) consider now a 
model equation 



de(r \t,r) 
dt 



= D T Ae{r Q \t,r). 



(2.13) 
Equa- 



which acco unts only for turbulent diffusion, 
tions ( |2 1CD and 02.1$ yield: 

SM,(t)/St=i9<|e(ro|t J r)|9- 1 i>rA|e(ro|t,r)|) 1< ,. (2.14) 



Subs titut ing distribution ( 2.11 ) w e esti mate Laplacian in 
Eq. 02.14) as -i/^d- Equations fl2.13p and 02.14) im p!y 
that 

dM q (t)/dt = -qD T M q (t)/e 2 cl . 
The solution of Eq. ( 2.15| ) reads: 

M q {t) = M g (0)exp[- 7di£ (gr)t], 



(2.15) 



7dif(<?) ~ qD T /i, 



cl - 



(2.16) 



where 7dif (?) denotes a contribution to the damping rate 
of M, q (t) caused by turbulent diffusion. 



2. Effect of particles inertia 

In this subsect ion we show that the term — O div W in 
the RHS of Eq. (|i]) can result in an exponential growth 
of M. q (t) oc exp[7i„(g) t], i.e., in the instability. We de- 
noted here the contribution to the growth rate of M. q (t), 
caused by the inertia of particles, by Ti n (?) ■ I n order 
to evaluate 7m(<?) we neglect now in Eq. 02. 8| ) both, the 
convective term in the LHS of this equation (i.e. the the 
turbulent velocity difference inside the cluster) and the 
molecular diffusion term. The resulting equation reads: 



de(r \t,r) 
dt 



= -Q(r Q \t,r) div W (r \t, r) . (2.17) 



The main contribution to the BL-velocity difference 
W(ro\t, r) in the RHS of this equation is due to the ed- 
dies with size £ c i, the characteristic size of the cluster. 
Denote by v c i the characteristic velocity of these eddies 
and by t v ~ l c \/v c \ the corresponding correlation time. 
In our qualitative analysis we neglect the r dependence 
of div W(ro|<, r) insid e the cluster and consider the di- 
vergence in Eq. 02.17 ) as a random process b(t) with a 
correlation time r„: 

divW(r \t,r) -> b(t) . (2.18) 

Together with the decomposition ( ^.11 ) this yields the 



following equation for the cluster amplitude A(t): 
8A{t) 



dt 



-A(t)b(t) . 



(2.19) 
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The solution of Eq. (2. IE) reads: 

A(t) = A exp[-I(i)] , I{t) = [ b{r)dT . (2.20) 

Jo 



Integral I(t) in Eq. (2.2C) can be rewritten as a sum of 
integrals /„ over small time intervals t„: 



*/t. n y 

I{t) = Y,In, Ut)= / b(r) 
n=l , ,s 



(n-l)r„ 



In our qualitative analysis integrals I n may be considered 
as independent random variables. Using the central limit 
theorem we estimate the total integral 



where (. . .) v denotes averaging over turbulent velocity 
ensemble, £ is a Gaussian random variable with zero 
mean and unit variance, N = t/r v . Now we calculate 



M q {t) = J ® q P(0 d( , P(C) = exp(-C 2 /2) • 

Therefore, M q (t) = J x exp(q 2 S 2 N/2), where 

Ji = (1/V2tt) J exp[-(C - qS^Nf/2] dC ~ 1 . 

Since the main contribution to the integral J\ arises from 
£ ~ qSyN, the parameter q cannot be large. In this 
approximation the q— moments 

M q (t)=M q {0)exp[y m {q)t] 

with 7i n (q) being the growth rate of the q-th moment due 
to particles inertia which is given by 

Tm(«) ~ l(Tv[divW(r \t,r)] 2 ) v q 2 . (2.21) 



3. Qualitative picture of the clustering instability 

In previous subsections we evaluated the contributions 
to the growt h rat e of M q (t) due to the turbulent diffusion 
7dif( <7), Eq . (2.16), and due to the particles inertia 7i n (<?)i 
Eq. ( [2.21 ). The total growth rate may be evaluated as a 
sum of these contributions: 

MS) = X g (0)exp( 7g i), (2.22) 

lq - 7dif(?) + 7in(<z) , 

lq qD T /e cl + ^(T v [dWW(r \t,r)} 2 ) v q 2 . 

Clearly, the instability is caused by a nonzero value of 
(r^(div W) 2 ), i.e., by a compressibility of the particle ve- 
locity field v(t, r). 

Compressibility of fluid velocity itself u(t, r) (including 
atmospheric turbulence) is often negligible, i.e., divu « 



0. However, due to the effect of particles inertia their 
velocity v(t, r) does not coincide with u(t,r) (see, e.g., 
p5| , p6| [2^ , pq|), and a degree of compressibility, a Vl of 
the field v(t, r), may be of the order of unity[|§ g ||. 
Parameter a v is defined as 



a v ee ([div«] 2 )/(|V 



X v\ 



(2.23) 



Note that parameter a v is independent of the scale of the 
turbulent velocity field. It characterizes a compressible 
part of the velocity field as a whole. The main contribu- 
tion to this parameter comes from the scales which are 
of the order of the Kolmogorov scale 77. 

For inertial particles divu ~ TpAP/p, where r p is the 
particle response time, 



t p = m p /6TTpva — 2p p a 2 /9pis , 



(2.24) 



where m p and p p are the mass and material density of 
particles, respectively. The fluid flow parameters are: 
pressure P, Reynolds number Tie = Lu^/v, the dissi- 
pative scale of turbulence 77 = LTZe~ 3 ^ 4 , the maximum 
scale of turbulent motions L and the turbulent velocity 
ut in the scale L. Now we can estimate a v as 



(p P /p) 2 (a/ V r = (a/a,) 4 



(2.25) 



(see 20 ), where a* is a characteristic radius of particles. 
For a > a, it is plausible to correct this estimate as 
follows: 

(2.26) 



a 



For water droplets in the atmosphere p p /p ~ 10 3 and 
a* ~ 77/3O. For the typical value of 77 ~ 1mm it yields 
a» ~ 30/im. On windy days when 77 decreases, the value 
of a* correspondingly becomes smaller. 

Then we estimate (r^[diw] 2 ) as 2a v /r v , because 
([iotv] 2 ) ~ 2t~ 2 . Assuming that the cluster size £ c \ is 
of the order of the inner scale of turbulence, 77, we have 
to identify t v with a turnover time of eddies in the inner 

scale 77, t v — > t v ee rj/Vrj = (L/u^Tie^ 1 ^ 2 . Thus, the 

growth rate 7 g of the g-th moment in Eq. ( 2.22) ) may be 
evaluated as 



lq - 7cl q(q - 9cr) , 7cl 



9cr - z^- ■ (2.27) 



Clearl y, the mome nts w ith q > q CI are unstable. Equa- 
tions ( 2.25 ) and ( 2.27 ) imply that it happens when 
a > a 9jCr where a q<CI = ai^/q 1 ^ 4 is the value of a at 
which q CT — q. The largest value of a ?jCr corresponds to 
the instability of the first moment, (|0|): ai, C r ~ 0.8a*, 
a 2 ,cr « 0.84ai,cr, a3,cr ~ 0.76ai,cr, a4cr ~ 0.71ai,cr, etc. 

Note that if (|0|) grows in time then almost all parti- 
cles can be accumulated inside the clusters (if we neglect 
a nonlinear saturation of such growth). We define this 
case as a strong clustering. On the other hand, if q CI > 1 
the first moment (|0|) does not grow, and the clusters 
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contain a small fraction of the total number of parti- 
cles. This does not mean that the instability of higher 
moments is not important. Thus, e.g., the rate of bi- 
nary particles collisions is proportional to the square of 
their number density (n 2 ) — (n) 2 + ( | G> | 2 ) - Therefore, 
the growth of the 2nd moment, (|0| 2 ), (which we de- 
fine as a weak clustering) results in that binary collisions 
occur mainly between particles inside the cluster. The 
latter can be important in coagulation of droplets in at- 
mospheric clouds whereby the collisions between droplets 
play a crucial role in a rain formation. The growth of 
the q-th moment, (|0| 9 ), results in that g-particles col- 
lisions occur mainly between particles inside the cluster. 
The growth of the negative moments of particles num- 
ber density (possibly associated with formation of voids 
and cellular structures) was discussed in |}0| (see also 

mil)- 

In the above qualitative analysis whereby we consid- 
ered only one-point correlation functions of the number 
density of particles, we missed an important effect of an 
effective drift velocity which decreases a growth rate of 
the clustering instability. For the one-point correlation 
functions of the number density of particles the effective 
drift velocity is zero for homogeneous and isotropic tur- 
bulence. However, in the equations for two-point and 
multi-point correlation functions of the number density 
of particles the effective drift velocity is not zero and as 
we will see in the next section it increases a threshold for 
the clustering instability. 

III. THE CLUSTERING INSTABILITY OF THE 
2ND MOMENT 

A. Basic equations 

In the previous section we estimated the growth rates 
of all moments (|6| 9 )- Here we present the results of a 
rigorous analysis of the evolution of the two-point 2nd 
moment 

$(*, R) = (6(t, r)0(t, r + R)) . (3.1) 

In this analysis we used stochastic calculus [e.g., Wiener 
path integral representation of the solution of the 
Cauchy problem for Eq. (^^), Feynman-Kac formula 
and Cameron-Martin-Girsanov theorem]. The compre- 
hensive description of this approach can be found in 

We showed that a finite correlation time of a turbulent 
velocity plays a crucial role for the clustering instability. 
Notably, an equation for the second moment $(t, H) of 
the number density of inertial particles comprises spatial 
derivatives of high orders due to a non-local nature of 
turbulent transport of inertial particles in a random ve- 
locity field with a finite correlation time (see Appendix A 
and [^0|). However, we found that equation for $(<, R) 
is a second-order partial differential equation at least for 
two models of a random velocity field: 



Model I. The random velocity with Gaussian statis- 
tics of the integrals f Q v(t' ,£)dt' and f*b(t',£)dt', see 
Appendix B. 

Model II. The Gaussian velocity field with a small yet 
finite correlation time, see Appendix C. 

In both models equation for <I>(t, R) has the same form: 

d<S>/dt = £${t,R), (3.2) 
L = B(R) + 2U(R) ■ V + D a/3 (R)V a Vp , 

but with different expressions for its coefficients. The 
meaning of the coefficients B(R), U(R) and D a p(R) is 
as follows: 

- Function B(R) is determined only by a compress- 
ibility of the velocity field and it causes generation of 
fluctuations of the number density of inertial particles. 

- The vector U(R) determines a scale-dependent drift 
velocity which describes a transfer of fluctuations of the 
number density of inertial particles over the spectrum. 
Note that U(R = 0) = whereas B(R = 0) ^ 0. For 
incompressible velocity field U(J?) = 0, B(R) = 0. 

- The scale-dependent tensor of turbulent diffusion 
D a p{R) is also affected by the compressibility. 

In very small scales this tensor is equal to the tensor 
of the molecular (Brownian) diffusion, while in the vicin- 
ity of the maximum scale of turbulent motions this ten- 
sor coincides with the usual tensor of turbulent diffusion. 
Tensor D a p(R) may be written as 

Dali{R) = 2DS aP + Dl p (R) , (3.3) 
Dl p (R) = Dl (O) - Dl (R) . 

In Appendix B we found that for Model I: 

B(R) » 2 / (6[0,€(r 1 |0)]6[r,€(r 2 |T)])dr, (3.4) 

JO 

rOO 

U(R) « -2/ («[0,£(ri|0)]&[T,£(r 2 |r)])dT, 
Jo 

/>oo 

Dl p {R) « 2 / («a[0,€(ri|0)]« /J [r,€(r 2 |r)])dT. 
Jo 

For the (^-correlated in time random Gaussian com- 
pressible velocity field the operator C is replaced by £o 
in the equation for the second moment $(i, R), where 

A, = B {R) + 2U (R) • V + D al3 (R)V a Vp , 
B (R) = V Q V /3 ^ Q/3 (i?) , (3.5) 
U , a (R) = V l3 D a p(R) 

(for details see [|o[ |2l[). In the ^-correlated in time ve- 
locity field the second moment $(f , R) can only decay in 
spite of the compressibility of the velocity field. The rea- 
son is that the differential operator £ = V ' Q V ' ijD a p{R) 
is adjoint to the operator C\ = b a p(R)\7 Q V p and their 
eigenvalues are equal. The damping rate for the equation 

d$/dt = £l$(t,R) (3.6) 
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has been found in Ref. j36| for a compressible isotropic 
homogeneous turbulence in a dissipative range: 



72 



(3 -ax) 2 



6r r) (l + ( TT)(l + 3a T ) 



Here ax is the degree of compressibility of the tensor 
D^pi^R). For the ^-correlated in time incompressible ve- 
locity field (cr T = 0) Eq. (|j§ was derived in Ref. [|§. 
Thus, for the Kraichnan model of turbulent advection 
(with a delta correlated in time velocity field) the clus- 
tering instability of the 2nd moment does not occur. 

A general form of the turbulent diffusion tensor in a 
dissipative range is given by 



D a JR) = {CiR 2 8 a t 



f C 2 R a Rp)/T ri , 

-I- cr T ) • 



(3.8) 



C\ =2(2 + ct t )/3(1- 
C 2 = 2(2a T - l)/3 (1 

The parameter ax is defined by analogy with Eq. ( 2.23 ) 
V £> T V V a Vf3Dlp(R) 



0~J> 



VxJDxxV V a VpDl, 0l {R)e aa , 1 ef 3p , 1 



(3.9) 



where e a g 7 is the fully antisymmetric unit tensor. Equa- 
tions ( 2. 23] ) and (3.9) imply that ox = <J V in the case 
of ^-correlated in time compressible velocity field. Equa- 
tions (^4) show that for a finite correlation time identi- 
ties ( |3.5| ) are violated and 

B(R) + B (R) , U(R)^U (R). 

For a random incompressible velocity field with a fi- 
nite correlation time the tensor of tu rbu lent diffusion 
D a0 {R) = T- 1 (C Q (r 1 )e /3 (r 2 )) [see Eq. @] and the de- 
gree of compressibility of this tensor is 



ox 



((v-O 2 ) 

((Vx£) 2 ) 



(3.10) 



where £(ri|t) is the Lagrangian displacement of a par- 
ticle trajectory which paths through point r% at t = 0. 
Remarkably that Taylor j3?|] obtained the coefficient of 
turbulent diffusion for the mean field in the form D^(R — 
0)=r- 1 (Uri\t)Uri\t)). 



B. Clustering instability in Model I 

Let us study the clustering instability for the model 
of the random velocity with Gaussian statistics of the 
integrals 



v(t',$dt'. 



b(t',£)dt', 



see Appendix B. In this model Eq. (3.2) in a non- 
dimensional form reads: 



<9$ 
~dt 



1 



m(r) lm(r) 
l/m(r) = (Ci +C 2 )r 2 +2/Sc, 



+ (u- c 2 y 



2$' 

r 



(3.11) 



where [/ = UR and Sc= z^/£> is the Schmidt number. 
For small incrtial particles advected by air flow Sc^S> 1. 
The non-dimensional variables in Eq. (3.11) are r = R/rj 



(3.7) and t = t/r^, B and U ar e mea sured in the units r, 



Consider a solution of Eq. (3.11) in two spatial regions. 

a. Molecular diffusion region of scales. In this region 
r <C Sc -1 / 2 , and all terms oc r 2 (with C\ , C% and U) may 
be neglected. Then the solution of Eq. (3.11) is given by 



where 



$(r) = (1 - ar 2 ) exp(7 2 i) , 



a = Sc(S- 72 t,)/12, B> l2 r n 



b. Turbulent diffusion region of scales. In this region 
Sc^ 1 / 2 < r « 1, the molecular d iffusi on term oc 1/Sc is 
negligible. Thus, the solution of ( 3.11 ) in this region is 



$(r) = A ir A exp( 72 i), 



where 



A = (C 1 -C 2 + 2U±iC 3 )/2(C 1 +C 2 ), 
C 2 = 4(B - 7 2 r 7) )(C 1 + C 2 ) ~ (C x -C 2 + 2U) 2 . 

Since the total number of particles in a closed volume is 
conserved: 



2 $(r) dr = 



This implies that C\ > 0, and therefore A is a complex 
number. Since the correlation function <I>(r) has a global 
maximum at r = 0, C\ > C 2 — 2f/. The latter condition 
for very small U yields <tt < 3. For r ^> 1 the solution 
for <I>(r) decays sharply with r. The growth rate 72 of 
the second moment of particles number density can be 
obtained by matching the correlation function <I>(r) and 
its first derivative $'(r) at the boundaries of the above 
regions, i.e., at the points r = Sc -1 / 2 and r = 1. The 
matching yields C 3 /2(C*i + C 2 ) » 2vr/lnSc. Thus, 



1 



72 



r n (l + 3cr T ) 
37r 2 (l + 3g T ) 2 ' 
(l + cr T )ln 2 Sc 



200cr u (cr T ~ O-u) (3-CTx) 2 



3(l + CTu) 6(l + cr T 

20(ct B - cru) 



7^(1 + CT B )(1 + f7u) 



(3.12) 



where we introduced parameters ctb and ou defined by 

B = 20<tb/(1 + <tb) , U- 20<ru/3(l + a v ) . (3.13) 

Note that the parameters ctb ~ cu ~ c„. For the 
(5-correlated in time random compressible velocity field 
ob = ctj = o~T = • Figure 1 shows the range of param- 
eters (a v , (Tx) for ctb = cru = o"« in the case of Sc = 10 3 
(curve c), Sc = 10 5 (curve b) and Sc — > oo (curve a). 
The dashed line cr„ = trx corresponds to the ^-correlated 
in time random compressible velocity field. This is a 
limiting line for the curve " a" . Figure 1 demonstrates 
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FIG. 1: The range of parameters (a v , or) for ob = ou = o\, 
in the case of Sc = 10 3 (curve c), Sc = 10 5 (curve b) and 
Sc — *• oo (curve a). The dashed line a v — <tt corresponds to 
the 5-correlated in time random compressible velocity field. 



that even a very small deviations from the ^-correlated 
in time random compressible velocity field results in the 
instability of the second moment of the number density 
of inertial particles. The minimum value of ctt required 
for the clustering instability is or ~ 0.26 and a corre- 
sponding value of a v w 0.12 (see Fig. 1). For smaller 
value of a v the clustering instability can occur, but it 
requires larger values of or- 

Notably, in Model II of a random velocity field (i.e., 
the Gaussian velocity field with a small yet finite correla- 
tion time) the clustering instability occurs when a v > 0.2 
(see Appendix C). Indeed, the growth rate 72 of the sec- 
ond moment of particles number density is determined 
by equation: 



B(a v )St 2 



(3 



6(1 



8(1 + 3cr„) / 7T 



B{a v ) 



3(1- 
12 ( b 2 



<Jv)(l 
2 



3<7„) 



<j v ) \lnSc / 
foai _ h \ 



where St= T ren /T?j is the Strouhal number, T = 72(1 
T rC n72) 2 , and 



«i 



2(19<t„ + 3) 
3(1 + a v ) 
1 

~27(l + <7„) 
850 / a v 



a 2 



2{ia v + 1) 
3(1 +a v ) ' 



-(12 - 1278cr„- 3067(7^) , 



9 Vl + cr^ 
1 

27(1 + a v ) 



(36 + 466^ + 2499CT 2 ) 




3 o 



FIG. 2: The range of parameters (Sc, a v ) for Sc = 10 3 (curve 
c), Sc = 10 5 (curve b) and Sc — > 00 (curve a). The line Sc = 
corresponds to the 5-correlated in time random compressible 
velocity field. 



Figure 2 shows the range of parameters (Sc, a v ) for Sc = 
10 3 (curve c), Sc = 10 5 (curve b) and Sc — > 00 (curve a). 

It is seen from FIG. |2| that for a v > 0.2 the second- 
order correlation function of the number density of iner- 
tial particles can grow in time exponentially (i.e., 72 > 0) 
even for very small Strouhal numbers. For example, in 
the vicinity of a v =3, the growth rate 72 of the cluster- 
ing instability of the second-order correlation function is 
given by 



72 = 4 x 10- 



n 



(3 - a v f 
240 



20 

y 



InSc 



(3.14) 



The sufficient condition for the exponential growth of the 
second moment of a number density of inertial particles 
is Sc > Sc*- 01 ), where the critical Schmidt number Sc^ 
is given by Sc*- 01 -* = Sc(72 = 0). The clustering instabil- 
ity occurs when the degree of compressibility of particles 
velocity a v > 0.2, i.e., for particl es an d droplets with the 
radius a* > 25.4 /mi. Equation ( [2.27 ) also yields a sim- 
ilar value <7 cr ~ 1/6 for the threshold of the inst abil ity 
of the 2nd moment (at q CT = 2). Note that Eq. (3.7) is 
written for Sc — > 00. 



IV. NONLINEAR EFFECTS 

The compressibility of the turbulent velocity field with 
a finite correlation time can cause the exponential growth 
of the moments of particles number density. This small- 
scale instability results in formation of strong inhomo- 
geneities (clusters) in particles spatial distributions. The 
linear analysis does not allow to determine a mechanism 
of saturation of the clustering instability. As can be 
seen from Eq. (3.12) molecular diffusion only depletes 



the growth rates of the clustering instability at the lin- 
ear stage (contrary to the instability discussed in Ref. 
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pO| ), The clustering instability is saturated by nonlinear 
effects. 

Now let us discuss a mechanism of the nonlinear satu- 
ration of the clustering instability using on the example of 
atmospheric turbulence with characteristic parameters: 
■q ~ 1mm, ~ (0.1 — 0.01)s. A momentum coupling of 
particles and turbulent fluid is essential when m p n c \ ~ p, 
i.e., the mass loading parameter <f> = m p n c \/ p is of the 
order of unity (see, e.g., [Q). This condition implies that 
the kinetic energy of fluid p(u 2 ) is of the order of the 
particles kinetic energy m p n c \(v 2 ) , where |u| ~ \v\. This 
yields: 



n c \ ~ a 3 (p/3p p ) 



(4.1) 



For water droplets p P /p ~ 10 3 . Thus, for a = a* ~ 
30/im we obtain n c i ~ 10 4 cm~ 3 and the total number 
of particles in the cluster of size r], N c i ~ ?7 3 ti c i ~ 10. 
This values may be considered as a lower estimate for the 
"two-way coupling" when the effect of fluid on particles 
has to be considered together with the feed-back effect 
of the particles on the carrier fluid. However, it was 
found in p8| that turbulence modification by particles is 
governed by the ratio of the particle energy and the total 
energy of the suspension (rather then the energy of the 
carrier fluid) and thus by parameter </>(l + </>) (rather then 
by <j) itself). Thus we expect that the two-way coupling 
can only mitigate but not stop the clustering instability. 

An actual mechanism of the nonlinear saturation of 
the clustering instability is "four way coupling" when 
the particle-particle interaction is also important. In 
this situation the particles collisions result in effective 
particle pressure which prevents further grows of con- 
centration. Particles collisions play essential role when 
during the life-time of a cluster the total number of colli- 
sions is of the order of number of particles in the cluster. 
The rate of collisions J ~ H c i/r^ can be estimated as 
J ~ 47ra 2 7ij!j|« re i|. The relative velocity v IC \ of collid- 
ing particles with different but comparable sizes can be 
estimated as \v te \\ ~ t p \(u ■ V)tt| ~ Tpui/r). Thus the 
collisions in clusters may be essential for 

n cl -a- 3 ( V /a)(p/3p p ), t s ~ a{Zap p /np) l '\ (4.2) 

where £ s is a mean separation of particles in the cluster. 
For the above parameters (a — 30/im) n c i ~ 3x 10 5 cm~ 3 , 
tg ~ 5 a w 150/mi and N c \ ~ 300. Note that the mean 
number density of droplets in clouds n is about 10 2 — 
10 3 cm -3 . Therefore the clustering instability of droplets 
in the clouds increases their concentrations in the clusters 
by the orders of magnitude. 

In all our analysis we have neglected the effect of sed- 
imentation of particles in gravity field which is essential 
for particles of the radius a > 100/xm. Taking t c \ ~ 77 
we assumed implicitly that r p < r^. This is valid (for 
the atmospheric conditions) if a < 60/im. Otherwise the 
cluster size can be estimated as l c \ ~ ^(Tp/r,,) 3 / 2 . 

Our estimates support the conjecture that the clus- 
tering instability serves as a preliminary stage for a co- 



agulation of water droplets in clouds leading to a rain 
formation. 



V. DISCUSSION 

In this study we investigated the clustering instability 
of the spatial distribution of inertial particles advected 
by a turbulent velocity field. The instability results in 
formation of clusters, i.e., small-scale inhomogeneities of 
aerosols and droplets. The clustering instability is caused 
by a combined effect of the particle inertia and finite cor- 
relation time of the velocity field. The finite correlation 
time of the turbulent velocity field causes the compress- 
ibility of the field of Lagrangian trajectories. The latter 
implies that the number of particles flowing into a small 
control volume in a Lagrangian frame does not equal to 
the number of particles flowing out from this control vol- 
ume during a correlation time. This can result in the 
depletion of turbulent diffusion. 

The role of the compressibility of the velocity field is 
as follows. Divergence of the velocity field of the inertial 
particles div v = r p AP/p. The inertia of particles results 
in that particles inside the turbulent eddies are carried 
out to the boundary regions between the eddies by iner- 
tial forces (i.e., regions with low vorticity and high strain 
rate). For a sm all molecular diffusivity divv oc —dn/dt 
[see Eq. Q)]. Therefore, dn/dt cx -T p AP/p. Thus 
there is accumulation of inertial particles (i.e., dn/dt > 0) 
in regions with AP < 0. Similarly, there is an out- 
flow of inertial particles from the regions with AP > 0. 
This mechanism acts in a wide range of scales of a tur- 
bulent fluid flow. Turbulent diffusion results in relax- 
ation of fluctuations of particles concentration in large 
scales. However, in small scales where turbulent diffu- 
sion is small, the relaxation of fluctuations of particle 
concentration is very weak. Therefore the fluctuations of 
particle concentration are localized in the small scales. 

This phenomenon is considered for the case when den- 
sity of fluid is much less than the material density p p of 
particles (p -C p p ). When p > p p the results coincide 
with those obtained for the case p -C p p except for the 



transformation r p — > /3*r p , where 



2 1 



P_ 

p P 



P P - P 
2p p + p 



For p > p p the value dn/dt oc — [3*t p AP / p. Thus there 
is accumulation of inertial particles (i.e., dn/dt > 0) in 
regions with the minimum pressure of a turbulent fluid 
since /3* < 0. In the case p > p p we used the equation of 
motion of particles in fluid flow which takes into account 
contributions due to the pressure gradient in the fluid 
surrounding the particle (caused by acceleration of the 
fluid) and the virtual ("added") mass of the particles 
relative to the ambient fluid (39| . 

The exponential growth of the second moment of a 
number density of inertial particles due to the small- 
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scale instability can be saturated by the nonlinear ef- 
fects (see Section IV). The excitation of the second mo- 
ment of a number density of particles requires two kinds 
of compressibilities: compressibility of the velocity field 
and compressibility of the field of Lagrangian trajecto- 
ries, which is caused by a finite correlation time of a ran- 
dom velocity field. Remarkably the compressibility of 
the field of Lagrangian trajectories determines the coeffi- 
cient of turbulent diffusion (i.e. the coefficient D a p near 
the second-order spatial derivative of the second moment 
of a number density of inertial particles in Eq. (3.2). 
The compressibility of the field of Lagrangian trajecto- 
ries causes depletion of turbulent diffusion in small scales 
even for a v — 0. On the other hand, the compressibility 
of the velocity field determines a coefficient B(r) near the 
second mom ent of a number density of inertial particles 
in Eq. (3.2). This term is responsible for the exponen- 
tial growth of the second moment of a number density of 
particles. 

Summary: 

• We showed that the physical reason for the cluster- 
ing instability in spatial distribution of particles in tur- 
bulent flows is a combined effect of the inertia of particles 
leading to a compressibility of the particle velocity field 
v(t,r) and a finite velocity correlation time. 

• The clustering instability can result in a strong clus- 
tering whereby a finite fraction of particles is accumu- 
lated in the clusters and a weak clustering when a finite 
fraction of particle collisions occurs in the clusters. 

• The crucial parameter for the clustering instability 
is a radius of the particles a. The instability criterion is 
a > a cr sa a* for which ((divu) 2 ) = (|rotv| 2 ). For the 
droplets in the atmosphere a* ~ 30/^m. The growth rate 
of the clustering instability 7 c i ~ r~ (a/a*) 4 , where t v 
is the turnover time in the viscous scales of turbulence. 

• We introduced a new concept of compressibility of 
the turbulent diffusion tensor caused by a finite corre- 
lation time of an incompressible velocity field. For this 
model of the velocity field, the field of Lagrangian tra- 
jectories is not divergence-free. 

• We suggested a mechanism of saturation of the clus- 
tering instability - particle collisions in the clusters. An 
evaluated nonlinear level of the saturation of the droplets 
number density in clouds exceeds by the orders of mag- 
nitude their mean number density. 
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APPENDIX A: Basic equations in the model with a 
random renewal time 



In this Appendix we derive Eq. (A2C) for the simul- 
taneous second-order correlation function $(i,r) which 
serves as a basis for further analysis in Appendixes [B| 
and |c| under some simplifying model assumptions about 
the statistics of the velocity field. 



Exact solution of dynamical equations for a 
given velocity field 

a. Simple case: no molecular diffusion 



Consider first Eq. ( p.l[ ) for the number density of par- 
ticles n(t, r) in the case D = 0: 



dn(t, r) 
dt 



V ■[n(t,r)v(t,r)} = 0, 



(Al) 



when all particles are transported only by advection. So- 
lution of Eq. ( |A1[ ) with the initial condition n(s,r) is 
given by 



n(t,r) = G(t,r)n[ S) &{t,r\s)], 



(A2) 



where £h(t,r\s) is the Lagrangian trajectory of the par- 
ticle which is located at coordinate r at time t. Here 
we label the particles at present moment of time t and 
consider a current time s < t as moments i n th e past. 
This differs from a usual approach, see Eqs. (2.4), when 
particles are labelled at the initial time to, and a cur- 
rent time t > t{). Therefore in the equations below it 
is more convenient to redefine Lagrangian displacement 
p L (t,r|s) -> p h (t,r\s) = -p h (t,r\s). Now Eqs. (|1|) can 
be written as 

p L (t,r\s) = J v[T,$ L (t,r\r)]dT, (A3) 
&(t,r\s) = r-p L (t,r\s) . (A4) 
The Green function is the functional of ^(t,r\s): 

G(t,r,s) = exp{- J &[r,£ L (i, r\r)] dr} , 

b(t,r) = V • v(t,r) , (A5) 
Introduce the shift operator 

exp[-p L • V] = 1 - p L • V + -[-pl ■ V] 2 - . . . (A6) 
which acts as follows: 

exp[-p L • V ]n(t, r) = n(t, r - p\) • (A7) 
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One can validate relation (A7) by Taylor series expan- 
sion of the function n(t,r — p^). Now Eq. (A2) can be 
rewritten as follows: 

n(i, r) = G(t, r, s) exp[-p L (t, r\s) ■ V }n(s, r) . (A8) 



b. Molecular diffusion as a Wiener process 



Consider now the full Eq. (2.1) with D ^ whereby 
particles are transported by both, fluid advection and 
molecular diffusion. It was found by Wiener (see, e.g., 
p3| ) that Brownian motion (molecular diffusion) can be 
described by the Wiener random process w(t) with the 
following properties: 



(w(t)) w = , (wi(t + r)wj(t)) w = rSij 



(A9) 



Here (. . .)„, denotes the mathematical expectation over 
the statistics of the Wiener process. Introduce the 
Wiener trajectory £w(*, T '|s) (which usually is called the 
Wiener path) and the Wiener displacement pw(i, r \ s ) as 
follows: 



£ w (t,r\s) = r-pw(t,r\s), 

pw(t,r\s) = J v[T,£ w (t, r|r)]dr + V2Dw 



(A10) 



Comparison of this formula with Eqs. ( |A3| ) shows that 
in the limit D — > 0, £w(i,f|s) — > and 
p w (t,r\s) -> p L (t,r\s). 

In Refs. J35| it was shown that solut ion of Eq. (|2~1 ) 
(with D 7^ 0) can be written as solution ( |A8|) of Eq. ( )A1 ) 
(with D = 0) by replacement p~\,{t,r\s) — > pw(t,r\s) 
and then averaging over the statistics of the Wiener pro- 
cesses (|A9| 



(All) 



i(t,r) = (G(t,r,s) exp[-pw(t,r\s) ■V]n(s,r)) l 



2. Two-step averaging over velocity statistics 



a. Model of a random velocity field 



Note that Eq. (All) is a solution of Eq. ( pTl| ) at a given 
realization of the random velocity field. Our next goal is 
to determine the simultaneous correlation functions 

n(t) = ({n(t,r))) v , (A12) 
$(t,r 2 -ri) = ((n{t,r 1 )n(t,r 2 ))) v -~n 2 (t) 7 

averaged over the stationary, space homogeneous statis- 
tics of turbulent velocity field, where ((. . .))„ denotes this 
averaging. Since the initial distribution n(to,r) is as- 
sumed to be homogeneous in space, n(t) is independent 
of spatial coordinate, and $(t, r 2 — ri) depends only on 
the difference r 2 — r±. 



In order to simplify the averaging procedure (A12) we 
consider a model of random velocity field which fully 
looses memory at some instants of renewal tj. For t\ 
and <2 inside a renewal interval [tj < ti,t% < t 3+ i] the 
velocity pair correlation function is defined as 

T afi (t 2 - t u r 2 - ri) = (v a (ti,ri)v (t 2 , r 2 )) v , (A13) 

where (. . .)„ denotes averaging over "intrinsic statistics" 
of the velocity field. In our model the velocity fields be- 
fore and after renewals are statistically independent. The 
interval between the renewal instants Tj may be the same 
or randomly distributed, say with the Poisson statistics. 
In the latter case the full averaging ((...))« may be con- 
sidered as a two-stage process. First one calculates (. . .)„ 
and then averages over the statistics of the renewal time 
T rcn , which is denoted as (. . . ) rcn ' 



(A14) 



For the Poisson statistics of t,- 



F^ih-tt^-n) = ((v a (t 1 ,r 1 )v (t 2 ,r 2 )))v 
= ^ Q/3 (t 2 -ti,r 2 -r 1 )exp(-|t 2 -t 1 |/f rcn ), (A15) 

where f ron is a mean renewal time. It would be useful to 
define the correlation time of the function as follows 



T V {R) 



(A16) 



Certainly this model of the random velocity field can- 
not be considered as universal. However, it reproduces 
important features of some flows (see, e.g., Ref. po|). 



b. Averaging procedure 

Our model involves three random processes: 

1. The Wiener random process which describes Brow- 
nian (molecular) diffusion. 

2. Poisson process for a random renewal time. 

3. The random velocity field between the renewals. 



Eq. ( AllD presents n(t,r) after the first step, i.e., it de- 
scribes the number density at a given realization of a 
velocity field. Using Eq. ( [AllD we obtain 

n(t, ri )n(t,r 2 ) = (G(n)G(r 2 ) exp^n) ■ V x 
+€'(r 2 ) ■ V 2 ]n(s, n)n(s, r 2 )) ww , (A17) 

where Vi = d/dr\ and V 2 = d/dr 2 and () ww de- 
notes averaging over two independent Wiener processes 
determining two Wiener paths. Hereafter for simplicity 
we use the following notations: G(r) = G(t,r,s) and 
e(r) = t'(t,r\s). 
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Now we average Eq. ( A17) over a random velocity field 
for a given realization of a Poisson process: 

$(*, r 2 - n) = (n(t, ri)n(t, r 2 )) v - [nf (A18) 
= «G(n)G(r 2 ) expK'(n) • V x + £'(r 2 ) ■ V a ])„„), 
X$(io,n - r 2 )} . 

Here the time to is the last renewal time before time t 
and t' = t — to is a random variable. Thus, averaging of 
the functions 

G(n)G(r 2 ) exp^n) • V x + £'(r 2 ) • V 2 ] , $(t , n - r 2 ) 

is decoupled into two time intervals because the first func- 
tion is determined by the velocity field after the renewal 
while the second function <5>(to,»"i — r 2 ) is determined 
by the velocity field before renewal. Now we take into 
account that for the Poisson process any inst ant c an be 
chosen as the initial instant. We average Eq. ( A18 ) over 
the random renewal time. The probability density pit!) 
for a random renewal time is given by 



P(*0 = T rcn exp(-t7 T ron) 



(A19) 



Thus the resulting averaged equation for "fully" averaged 
corre latio n function $(t,i£) = i?)) rcn , defined by 

Eq. ( A12| ), assumes the following form: 



. P(T,R)$(t- T,R)cxp(-T/T ven )dT 
10 



■exp(-t/f ren )P(t,J2)$o(-R) 



(A20) 



The first term in Eq. ( A20 ) describes the case when there 
is at least one renewal of the velocity field during the 
time t (i.e., the Poisson event), whereas the second term 
describes the case when there is no renewal during the 
time t. Here $ (-R) =$(t = 0,R) and 



P(t,R) = «G(n)G(r 2 ) 

x exp[£'(j"i) • Vi + £'(r 2 ) • V 2 
exp((g( ri ) + g{r 2 ) + £'(n) • V x + £'(r 2 ) • V 2 ) 



where G(r) = exp[g(r)]. Equation ( A20 ) is simplified in 
Appendixes B and C under the additional assumptions 
about the velocity field statistics. 



APPENDIX B: VELOCITY FIELD WITH 
GAUSSIAN LAGRANGIAN TRAJECTORIES 

Consider a model of a random velocity field where La- 
grangian trajectories, i.e., the integrals J v((i, £) d(i and 
J b{(i, £) d(i have Gaussian statistics. Using an identity 



(exp(ary))„ = exp(ia 2 ) in Eq. ( A21) we obtain 



where 

t = B(R)+2U a {R)V a +D aP {R)V a Vf3, 

(iB{R) = ((g(r 1 )g(r 2 )) ww ) v , (B2) 

(iU a (R) = -{{£' a (ri)g(r 2 )) ww ) v , 

D a0 (R) = D a0 (O) - D a0 (R) , 

liD a p(R) = ((C(ri)^(r 2 )) ww ) v . 

Here rj is a Gaussian random variable with zero mean 
value and unit variance and ((G(r)) w ) v = 1. The latter 
yields 

((g) w )v = -^((g 2 )zu)v , g = ((g) w ) v + g , 

wher e (( g) w ) v = 0. When correlation time t v (R), 
Eq. ( Alq ), is much less then the current time t and f rcn , 
these correlation functions are given by 



B(R) = 2j ((b[0^( ri )}b[ f /^(r 2 )}) ww ) v d^, (B3) 



00 

U a {R) = -2j ({v a %£{ri)]b\p! ,£{r 2 )]) ww ) v dii' , 



00 

D af3 (R) = 2j {(v a [O^( ri )}v [(i'^(r 2 )}) ww ) v d(i' , 


where we used an identity 

(( / a a ((i' ,Ti) dp,' / cp{(i" ,r 2 ) d(i") 



2(i / ((a a (0,ri)cf3((i' ,r 2 )) w ) v d(i' 



(A21) auows t° rewrite Eq. ( A20| ) as 



$(t,R) = J- 



exp(fiLi) d(i 



<f>(t,R) 



where 



+ exp(t£ 1 )<f>(t,R), 



dt r ron 



(B4) 



(B5) 



To derive Eq. (B4) we used the following identity 

S(f - (i, R) = exp f-A*^) *(*, R) > (B6) 
which follows from the Taylor expansion 



P{(i, R) = exp[^£] , 



(Bl) 



/(^) = £(4) 

m— 1 v / 
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In particular, 



$o(H) = *(* -t,R) = exp ( -t— ) *(t, R) 



d 



Evaluating the integral in Eq. (B4) we obtain 

exp(tA) - l] (ti + f-^j R) = 

Here we used the commutativity relation 
£1 exp(i£i) = exp(tCi) L\ . 

Thus, finally 



<9$ 



B(R) + 2V(R) ■ V + D al3 {R)V a V \ i?) , 

(B9) 



(B8) 



Note that in the limit f rcn — > oo, Eq. (BE) describes the 
evolution of <!>(£, R) in the model of the random velocity 
field without renewals. 



where M = l + T rcn (<9/<9£) and we neglected the last term 
in Eq. ( A20| ) for small r. Expanding the function P(t, R) 
in Taylor series in the vicinity of r = we obtain 



oo 
k=0 



d k P(r, R) 



Q T k 



£, r {k+i) _ jj) = o, 

t— o > 



(C2) 



where we used that 



/ 



T k exp ( __m ) dr = k\ fi+ 1 M-( fe+1 ) 



Neglecting the terms ~ 0(r r 5 cn ) in Eq. (|C2]) we obtain 



APPENDIX C: GAUSSIAN VELOCITY FIELD 
WITH A SMALL YET FINITE CORRELATION 
TIME 

Here we consider a random Gaussian velocity field with 
a small f rcn . Using Eq. (B6) we rewrite Eq. ( |A20| ) in the 
form 



M 



a d*(t,R) 



Of 



d 2 P(r,R) 

' d*P(T,R) \ 

<9r 4 



(C3) 



T = 0~ 



*(*,«), 



■J— I P(T,R)exp(-^—M )dr-l }<P(t,R) = 0, 

7"rcn J \ ^"rcn 







(Cl) 



where we used that the expansion of the operator P{t, R) 
into Taylor series (for small r) for a random Gaussian ve- 
locity field has only even powers of r. Thus, the equation 
for the correlation function $(i, R) is given by 



M 2 mt dt R) = [B{R)+2U{R)-V + b a p{R)V a V ]$, (C4) 

where 

D a p(R) = ^((^G( ri )G(r 2 )) ww ) v , (C5) 

U a (R) = -J-{{g(r 2 )e a (ri)) ww ) v + ^{{g(ri)g(r 2 )Z a ) ww ) v , (C6) 

5(H) = ^-((g(r 1 ) 9 (r 1 )) ww ) v , (C7) 

Tren 

Here for the homogeneous turbulent velocity field plj : 

i = ? V = d/dR, G = G + g, ((g) ww ) v =Q, G = ((G) ww ) v = 1 . 

Using the expansion of £(T ren ,r) and g[T re n, into Taylor series of a small time f rcn after the lengthly algebra we 
obtain 

D a p{R) = 2£>J a/3 + 2f ron [/ Q/3 (i?) + St 2 Q a0 {R)} , (C8) 
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Q a p{R) = 3[(V l/ /^)(V M / Q ,) - ftutVuVpfaf}] + 2AA a Ap + 12(i4 M V^/ a/ j - /a/jV^) - 20/ Q/1 V M A /3 , (C9) 
U a (R) = -2f I , n {A Q -St 2 [(V^ A1 )(V M / cw ) + 10^V Al A ct + 12^ Q V M A Al ]} , (CIO) 
B(R) = -2f rcn {V AI A M + St 2 [(V^)(V M ^)-6(V M A M ) 2 ]} , (Cll) 

= ^pfap, fa/3 = fap(0) - fap(R) , fap(R) = (v a (ri)Vp{r 2 ))v , 



and St= T Icn /T v is the Strouhal number. In these cal- 
culations we neglected small terms ~ 0(St 2 i? 3 V 3 ). Our 
analysis showed that the neglected small terms do not 
affect the g rowt h rate of the clustering instability. In 
Eqs. (|C9|)-( CI]) ) we assumed that the correlation func- 
tion f a fj for homogeneous, isotropic and compressible ve- 
locity field is given by 



,2 r 



fafl(R) 



RF' 

{F + F c )8 aP + —P Q/ 3 + RF' c R a p 



(C12) 

(see ]36|), and in scales < R <C 1 incompressible F(R) 
and compressible F C (R) components of the random ve- 
locity field are given by 

F(R) = (1 - i? 2 )/(l + o w ) , F C {R) = a v F(R) , 

in scales R > 1 the functions F — F c = 0. Here R 
is measured in the units of 77, P a p{R) — &ap — Rap, 
R a0 = R a Rp/R 2 and F' = dF/dR. Turbulent diffusion 
tensor D a p (R) is determined by the field of Lagrangian 
trajectories £ [see Eq. (|C5|)1 . Due to a finite correla- 
tion time of a random velocity the field of Lagrangian 
trajectories £ is compressible even if the velocity field is 
incompressible (a v = 0). Indeed, for a v — we obtain 



«(v-0 2 



Using Eqs. Q-(|Clf) we calculate the functions 
D a p(R), U a (R) and B(R) : 



D a p(R) = [2D + R 2 {a 3 + St 2 b 6 )}S a p 
+R 2 (a 4 + St 2 b 4 )R aP , 
U a (R) = -R a (a 5 + St 2 b 5 ), 

B = 0,6 + St 2 & 2 , 



(C13) 
(C14) 
(C15) 



where b 2 



«5 



«4 



b 6 = 



b 4 



%h and 



20er 1 , 
~3(l + <7„) 
4o„-2 
3(1 + o„) ' 
12 + 872o„ 



- -^1 

3 

h = - 

- 433o 2 



a-3 

2350 
27 V 1 



2a v 



3(l + o-«) 
2 



27(1 + a v ) 2 
2(12- 203o„ + 1033og) 
27(1 + <7„) 2 



We will show here that the combined effect of particles in- 
ertia (a v ^ 0) and finite correlation time (St 7^ 0) results 



in the excitation of the clustering instability whereby un- 
der certain conditions there is a self-excitation of the sec- 
ond moment of a number density of inertial particles. 
This instability causes formation of small-scale inhomo- 
geneities in spatial distribution of inertial particles. 

The equation for the second-order correlation function 
for the number density of inertial particles reads 



M 



2 d$(t,R) _ $' 



at 



m(R) 



\(R)& + B<$> 



(C16) 



see Eqs. (C15)], where the time t is measured in units of 

2(1 + X 2 ) 



t v , and 



A 

P 
01 

61 

&3 



dR' 



<9 2 $ 



2[2 + X 2 (l + 2C)] 



RSc 



1 

m 

C = 



a 1 



Sc 

+- St 2 6i 



4/i 



a 2 + St 2 6 3 



X(R) = y/ScfiR, R=\r 2 -ri\ 



2(19o„ + 3) 
3(1 + a v ) 
1 

"27(1- 



0,2 



2(3gg + 1) 
3(l + o,) ' 

-2(12- 1278o„- 3067o- 2 ) , 



27(1 + a v y 



-(36 + 4660^ + 2499o 2 ) 



In order to obtain a solution of Eq. ( C16| ) we use a sep- 
aration of variables, i.e., we seek for a solution in the 
following form: 

$(*, R) = &(R) exp( 72 i) , 

whereby 72 is a free parameter which is determined using 
the boundary conditions 



$(i? = 0) = l, $(#->• 00) 



Q 



Here 72 is measured in units of l/t v . Since the func- 
tion $(i, R) is the two-point correlation function, it has 
a global maximum at R — and therefore it satisfies the 
conditions: 

<£'(# = 0) = 0, $"(i? = 0) < 0, 

$(i? = 0) > \$(R > 0)| . 



Then Eq. ( |Ciq ) yields 
1 



r$(i?) 



m(i?) 



4" + \(R)& + B$ , 



(C17) 
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where T = 72 (1 + r rcn -f 2 ) 2 
solution for < R < 1 : 



Equation (C17) has an exact 



$(X) = S(X)X(1 + X 2 )^ 2 , 

S(X) = Be{A 1 P^{iX) + A 2 Q!i{iX)} 1 



(C18) 



and P£(Z) and Q*(Z) are the Legendre functions with 
imaginary argument 

B-T 



z = ix, v = c-^, ( = -±±Vc 2 -k, K = , ); 



Solution of Eq. flClq ) can be analyzed using asymp- 
totics of the exact solution (C18). This asymptotic analy- 
sis is based on the separation of sca les (s ee, e.g., J34], [36f ). 
In particular, the solution of Eq. (C16) has different re- 
gions where the form of the functions m(R) and A(i?) are 
different. The functions $(i?) and in these differ- 

ent regions are matched at their boundaries in order to 
obtain continuous solution for the correlation function. 
Note that the most important part of the solution is lo- 
calized in small scales (i.e., R <§; 1). Using the asymptotic 
analysis of the exact solution for X ^> 1 allowed us to ob- 
tain the necessary conditions of a small-scale instability 
of the second moment of a number density of inertial par- 
ticles. The results obtained by this asymptotic analysis 
are presented be low. 

The solution (CIS) has the following asymptotics: for 
X <C 1 (i.e., in the scales < R <C 1/vSc) the solution 
for the second moment $ is given by 



Hx) = {l-(n/6)[X 2 + 0(X i )}} 



(C19) 



For X > 1 (i.e., in the scales 1/vSc -C R < 1) the 
function $ is given by 



= Re{AX- c±vc ' 2 - K }. 



(C20) 



When C 2 — k < the second-order correlation function 
for a number density of inertial particles $ is given by 



$(i?) = A Z R~ C cos{v! In R + tp) 



where C > and (p is the argument of the complex con- 
stant A. For R > 1 the second-order correlation function 
for the number density of inertial particles is given by 



$(i?) = (A4/R) exp(-R VST/2) , 



(C21) 



where T > 0. Since the total number of particles in 
a closed volume is conserved, i.e., particles can only be 
redistributed in the volume. 



f 

Jo 



R 2 <f>(R) dR = $(fc = 0) = 



The latter yields <p = -it/2 fo r InS c > 1 and r« 1. 
When C 2 — k > 0, t he sol ution (C2C) cannot be matched 
with solutions (jCl9D and ( |C2l[) . Thus, the condition C 2 - 
k < is the necessary condition for the existence of the 
solution for the correlation function. The condition C > 
provides the existence of the global maximum of the 
correlation function at R = 0. 

Matching functions $ and $' at the boundaries of the 
above-mentioned regions yields coefficients Aj» an d T. In 
particular, the eigenvalue T is given by Eq. (|3.14 ). 
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